Primordial Non-Gaussianity in the Cosmic Microwave Background 



Amit P.S. YadavlB and Benjamin D. Wandelt 

^ Institute for Advanced Study, School of Natural Sciences, Einstein Drive, Princeton, N J 08540, USA 
^ UPMC Univ Paris 06, Institut d'Astrophysique de Paris, 98 bis, boulevard Arago 75014 Paris, France and 
Department of Astronomy and Physics, University of Illinois at Urbana-Champaign, Urbana IL 61801, USA 

In the last few decades, advances in observational cosmology have given us a standard model of cosmology. 
We know the content of the universe to within a few percent. With more ambitious experiments on the way, 
we hope to move beyond the knowledge of what the universe is made of, to why the universe is the way it 
is. In this review paper we focus on primordial non-Gaussianity as a probe of the physics of the dynamics of 
the universe at the very earliest moments. We discuss 1) theoretical predictions from inflationary models and 
their observational consequences in the cosmic microwave background (CMB) anisotropies; 2) CMB-based 
estimators for constraining primordial non-Gaussianity with an emphasis on bispectrum templates; 3) current 
«j ' constraints on non-Gaussianity and what we can hope to achieve in the near future; and 4) non-primordial 

I sources of non-Gaussianities in the CMB such as bispectrum due to second order effects, three way cross- 

correlation between primary-lensing-secondary CMB, and possible instrumental effects. 

I. MOTIVATION 

o 

' In the last few decades the advances in observational cosmology have led the field to its "golden age." Cosmologists are 
beginning to nail down the basic cosmological parameters. We now know that we live in a Universe which is 13.7 ± 0.1 Gyr old 
1^ and is spatially flat to about 1%, and is made of 4.6 ± 0.1% baryons, 22.8 ± 1.3% dark matter, and remaining 72.6 ± 1.5% in the 
form of dark energy. Although we know the constituents to high accuracy, we still do not completely understand the physics of 
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Q ■ the beginning, the nature of dark energy and dark matter Many upcoming CMB experiments complimented with observational 
|H ' campaign to map 3D structure of the Universe and new particle physics constraints from the Large Hadron Collider will enable 
^ . us to move beyond the knowledge of what the universe is made of, to why the universe is the way it is. In this paper we focus on 
I I ' learning about the physics responsible for the initial conditions for the universe. 

Inflation 111-3] is perhaps one of the most promising paradigms for the early universe, which apart from solving some of 
CO ■ the problems of the Big Bang model like the flatness and horizon problem, also gives a mechanism for producing the seed 
' perturbations for structure formation IH-l^, and other testable predictions 

Most observational probes based on 2-point statistics like CMB power spectrum still allow vast number of inflationary 
models. Moreover, the alternatives to inflation such as cyclic models are also compatible with the data. Characterizing the 
> non-Gaussianity in the primordial perturbations has emerged as powerful probe of the early universe. The amplitude of non- 
Gaussianity is described in terms of dimensionless non-linearity parameter /nl (defined in SecHIIli. Different models of inflation 
predict different amounts of /nl, starting from 0(1) to /nl ^ 100, above which values have been excluded by the WMAP data 
already. Non-Gaussianity from the simplest inflation models that are based on a slowly rolling scalar field is very small lfloljl5ll : 
however, a very large class of more general models with, e.g., multiple scalar fields, features in inflaton potential, non-adiabatic 
fluctuations, non-canonical kinetic terms, deviations from Bunch-Davies vacuum, among others lfl6l for a review and references 
therein] generates substantially higher amounts of non-Gaussianity. 

The measurement of the bispectrum of the CMB anisotropies is one of the most promising and "clean" way of constraining 
i ' /nl- Many efficient methods for evaluating bispectrum of CMB temperature anisotropies exist [fl7l - l2lll . So far, the bis- 
pectrum tests of non-Gaussianity have not detected any significant /nl in temperature fluctuations mapped by COBE [|22ll 
and WMAP lfT8ll23l - 58ll . On the other hand, some authors have claimed non-Gaussian signatures in the WMAP temperature 
data ll29l - l33ll . These signatures cannot be characterized by /nl and are consistent with non-detection of /nl- 

Currently the constraints on the /nl come from temperature anisotropy data alone. By also having the polarization infor- 
mation in the cosmic microwave background, one can improve sensitivi ty to p rimordial fluctuations ifsil Issll . Although the 
experiments have already started characterizing polarization anisotropies II36I439I1 . the errors are large in comparison to temper- 
ature anisotropy. The upcoming experiments such as Planck will characterize polarization anisotropy to high accuracy. 

The organization of the paper is as following: In Section we review the inflationary cosmology focusing on how the 
microscopic quantum fluctuations during inflation gets converted into macroscopic sees perturbations for structure formation, 
and as CMB anisotropies. In Section|III]we discuss theoretical predictions for non-Gaussianity from the inflationary cosmology. 
In Section |IV] we show how the primordial non-Gaussianity is connected to the CMB bispectrum, and describe/review CMB 
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bispectrum based estimators to constrain primordial non-Gaussianity (/nl)- In Section |V] we discuss the current constraints 
on /nl by CMB bispectrum and what we can hope to achieve in near future. We also discuss non-primordial sources of non- 
Gaussianity which contaminate primordial bispectrum signal. In section |VT] we discuss other methods for constraining /nl 
besides CMB bispectrum. Finally in Section lVTIl we summarize with concluding remarks. 



II. INTRODUCTION: THE EARLY UNIVERSE 

One of the most promising paradigms of the early universe is inflation d 0. 01, which apart from solving the flatness, 
homogeneity and isotropy problem, also gives a mechanism for producing the seed perturbations for structure formation, and 
other testable predictions ' (for a recent review of inflationary cosmology see [43]). During inflation, the universe goes through 
an exponentially expanding phase. From the Friedman equation, the condition for the accelerated expansion is 

p + 3p<0. (1) 

For both matter and radiation this condition is not satisfied. But it turns out that for a scalar field, the above condition can be 
achieved. For a spatially homogeneous scalar field, (j>, moving in a potential, V{(f)), the energy density is given by 

p^ = ^(d(t>/dt)^ +Vi(P), (2) 

p^^lid(f>/dt)^-V{4>). (3) 

Hence the condition for accelerated expansion of the universe dominated with scalar field (p is 

{dcj)/dt)^ < y (</)). (4) 

Physically this condition corresponds to situations where kinetic energy of the field is much smaller than its potential energy. 
This condition is referred to slowly-rolling of the scalar field. During such slow-roll, the Hubble parameter, H{t) = d In a/dt, is 
nearly constant in time, and the expansion scale factor, a{t) is given by 

a{t) = a(to) exp ( [ H{t')dt'\ w a(io) cxp [H{t) {t - to)] . (5) 



and the pressure is given by 



This exponential expansion drives the observable universe spatially flat, homogeneous and isotropic. 

A toy model is shown in Fig.[T] In the slow-roll phase, rolls down on V{(t)) slowly, satisfying Eq. (01) and hence driving the 
universe to expand exponentially. Near the minima of the potential, ((> oscillates rapidly and inflation ends. After inflation ends, 
interactions of (j) with other particles lead (f) to decay with a decay rate of F^, producing particles and radiation. This is called a 
reheating phase of the universe, as cj) converts its energy density into heat by the particle production. 

Not only inflation solves the flatness, homogeneity and isotropy problem, it also gives a mechanism for generating seed 
perturbations. During inflation the quantum fluctuation in the field (p are exponentially stretched due to the rapid expansion 
phase. The proper wavelength of the fluctuations are stretched out of the Hubble-horizon scale to that time, H~^. Once outside 
the horizon, the characteristic rm.s. amplitude of these fluctuations is |'?!>|;,jns ~ H/{2tt). These fluctuations do not change in 
time while outside the horizon. After inflation, and reheating, the standard hot-big scenario starts. As the universe decelerates, 
at some point the fluctuations re-enter the Hubble horizon, seeding matter and radiation fluctuations in the universe. Figure |2] 
summarizes the evolution of characteristic length scales. 



Although inflation is the most popular theory for the early universe, other mechanisms, for example, ekpyrotic models (40|| and cyclic models l4lll43l have 
been proposed for generating nearly scale invariant Gaussian perturbations, while retaining homogeneity and flatness. In the cyclic universe, there is no 
beginning of time, and our expansion of the universe is one out of the infinite number of such cycles. Each cycle consists of the following phases: (1) A hot 
big bang phase, during which a structure formation takes place. (2) An accelerated expansion phase which dilutes the matter and radiation energy density. 
Since observations suggest that our universe is going through an accelerated expansion phase, in the cyclic model interpretation, we are presently going 
through this phase. (3) A decelerating phase, which makes the universe flat, and generates nearly Gaussian and scale invaiiant density perturbations. (4) A big 
crunch/bang transition phase during which matter and radiation is created. Although the mechanism is different, the outcome of phase (3) of the cyclic model 
is in some sense analogous to a slow-roll expansion phase of inflation; and phase (4) will coiTespond with the reheating phase in the inflationary scenario. As 
we will discuss in the next section these two scenarios can be distinguished by their different predictions about the gravitational waves, and non-Gaussianity. 
Cyclic models predict negligible contribution of gravitational waves while inflationary models can produce large gravitational wave contribution, which can be 
detected by next generation experiments. Second, cyclic models produce much larger non-Gaussianity (of local type) in comparison to the standard slow-roll 
inflationary scenario. 
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Primordial Perturbations 



We use linearly perturbed conformal Friedmann Lematre Robertson Walker (FLRW) metric of the form, 

ds^ = a^ir) {-(1 + 2AQ)dT^ - 2BQdTdx' + [(1 + 2Hi^Q) + 2HTQij] dx'dx^} 



(6) 



where all the metric perturbations. A, B, H]^, and H^y, are ^ 1, and functions of conformal time r. 
dependence of the perturbations is described by the scalar harmonic eigenfunctions, Q, Qi, and Qij, 



-k'^Q, Q, = -k ^Q^i, and = k "^Q^ij + ^Si^Q 
Lets consider two new perturbation variables ||8I144 



Note that Qij is traceless: S^^ Qi 



0. 



The spatial coordinate 
that satisfy S^^Q^ij = 



and 



u = 6(t) -7^, 

ah 



aH aH ^ ^ 
C = —u = TZ —0(p, 



(7) 



(8) 



which are Gauge invariant. Here TZ = Hi, + ^Ht, is perturbations in the intrinsic spatial curvature. While u reduces to 6(j) 
in the spatially flat gauge (TZ = 0), or to ~{(j)/aH)TZ in the comoving gauge {5(t> = 0), its value is invariant under any gauge 
transformation. Similarly C,, which reduces to TZ in the comoving gauge, and to — {aH / (f))5(f) in the spatially flat gauge, is also 
gauge invariant. The perturbation variable C, helps the perturbation analysis not only because of being gauge invariant, but also 
because it is conserved on super-horizon scales throughout the cosmic evolution. 

The quantum fluctuations generate the gauge-invariant perturbation, u, that reduces to either 5(t) or {(f)/aH)TZ depending on 
which gauge we use, either the spatially flat gauge or the comoving gauge. Hence, Scfinat and {(f) / aH)TZcom are equivalent to each 
other at linear order. The benefit of using u is that it relates these two variables unambiguously, simplifying the transformation 
between (J^dat and T^com- 

The solution for ( is vahd throughout the cosmic history regardless of whether a scalar field, radiation, or matter dominates 
the universe; thus, once created and leaving the Hubble horizon during inflation, ( remains constant in time throughout the 
subsequent cosmic evolution until reentering the horizon. The amplitude of C is fixed by the quantum-fluctuation amplitude in u 
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2'!l{d(j)/dt) 



(9) 



This is the spectrum of C, on super-horizon scales. 
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FIG. 2: Evolution of comoving horizon and generation of perturbations in the inflationary universe. Figure from Ref. 143 



From Primordial Perturbations to CMB Anisotropies 



The metric perturbations perturb CMB, producing the CMB anisotropy on the sky. Among the metric perturbation variables, 
the curvature perturbations play a central role in producing the CMB anisotropy. 

As we have shown in the previous subsection, the gauge-invariant perturbation, does not change in time on super-horizon 
scales throughout the cosmic evolution regardless of whether a scalar field, radiation, or matter dominates the universe. The 
intrinsic spatial curvature perturbation, TZ, however, does change when equation of state of the universe, w = p/p, changes. 
Since C remains constant, it is useful to write the evolution of TZ in terms of C and w; however, TZ is not gauge invariant itself, 
but C is gauge invariant, so that the relation between TZ and ^ may look misleading. In 1980, Bardeen 14511 introduced another 
gauge-invariant variable, $ (or <i>H in the original notation), which reduces to TZ in the zero-shear gauge, or the Newtonian 
gauge, in which B = = Ht- $ is given by 



aH 



-B 




(10) 



Here, the terms in the parenthesis represent the shear, or the anisotropic expansion rate, of the t = constant hypersurfaces. 
While $ represents the curvature perturbations in the zero-shear gauge, it also represents the shear in the spatially flat gauge in 
which TZ = Q. Using we may write C, as 



aH 

C = TZ ^d<j) = $ 




(11) 



where the terms in the parenthesis represent the gauge-invariant fluid velocity. 

We use $ in rest of the paper because it gives the closest analogy to the Newtonian potential, which we have some intuition of. 
$ reduces to TZ in the zero-shear gauge (or the Newtonian gauge) in which the metric (Eq.©) becomes just like the Newtonian 
limit of the general relativity. 

The gauge-invariant velocity term, v — fc^^iJT, differentiates ( from $. Since this velocity term depends on the equation of 
state of the universe, w ^ p/p, the velocity and $ change as w changes, while ( is independent of w. The evolution of $ on 
super-horizon scales in cosmological linear perturbation theory gives the following ll46ll . 



(12) 
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for adiabatic fluctuations, and hence $ = in the radiation era (w = 1/3), and $ = |C in the matter era (w = 0). $ then 
perturbs CMB through the so-cafled (static) Sachs-Wolfe effect 147i| 

AT 1 + w 



T 5 + 3m; 



C- (13) 



At the decoupHng epoch, the universe has already been in the matter era in which w = 0, so that we observe adiabatic 
temperature fluctuations of AT/T = — ^'I' = ^^C' the CMB fluctuation spectrum of the Sachs-Wolfe effect, Ag^(fc), is 



Aiw(fc) = Tl^Uk) = ^Alik) (14) 



iA|(fc) = - 
9 *^ ^ 25 

By projecting the 3-dimension CMB fluctuation spectrum, A|-yy (fc), on the sky, we obtain the angular power spectrum^, Ci 14 

.sw_.„/ rffc A2 n.^,2,u^ ^ ^swr[(9-n,)/2]r[Z + (n,-l)/2] 



Jo k swv 7J/ L V dcc;j 2 r[(n^ + 3)/2]r[/ + (5-n,)/2]' ^ ^ 

where tq and Tdcc denote the conformal time at the present epoch and at the decoupling epoch, respectively, and Ug = I + 
[din A^(fc) /din fc] is a spectral index which is conventionally used in the literature. 

On small angular scales (£ > 10), the Sachs- Wolfe approximation breaks down, and the acoustic physics in the photon-baryon 
fluid system modifies the primordial radiation spectrum ll49tl . To calculate the anisotropics at all the scales, one has to solve the 
Boltzmann photon transfer equation together with the Einstein equations. These equations can be solved numerically with the 
Boltzmann code such as CMBFAST IdOII . The CMB power spectrum then can be written as 



= 471 / — A|(%|,(fc). (16) 
Jo k 

Here gTi{k) is called the radiation transfer function, and it contains all the physics which modifies the primordial power 
spectrum A$ to generate CMB power spectrum C^. For the adiabatic initial conditions, in the Sachs- Wolfe limit, gTi{k) = 
Ji [^■('''o ^ ''dec)]- Often in the literature power spectrum, P$(fc), is used instead of A|(A;). The two are related as A|,(fc) = 
(27r^)~^fc'^P$(fc). A|(fc) is called the dimensionless power spectrum. 

If $ were exactly Gaussian, all the statistical properties of $ would be encoded in the two-point function or in C( in the 
spherical harmonic space. Since $ is directly related to C through Eq. ( fT2] l. all the information of C is also in-coded in Cg. 
Although C which is related to a Gaussian variable, u, through C = —{aH/<j>)u, in the linear order ( also obeys Gaussian 
statistics; however the non-linear relation between ( and u makes ( (and hence $ and CMB anisotropics) slightly non-Gaussian. 
The non-linear relation between C and <i> is not the only source of non-Gaussianity in the CMB anisotropics. For example, at the 
second order, the relationship between $ and AT/T is also non-linear. 



Probes of the Cosmological Initial Conditions 

The main predictions of a canonical inflation model are: 

• spatial flatness of the observable universe, 

• homogeneity and isotropy on large angular scales of the observable universe, 

• seed scalar and tensor perturbation with primordial density perturbations being 

(a) nearly scale invariant, 

(b) nearly adiabatic, and 

(c) very close to Gaussian. 

At the time of writing, these predictions are consistent with all current observations. This represents a major success for the 
inflationary paradigm. On the other hand, the inflationary paradigm can be realized by a large 'zoo' ^ of models. In addition. 



2 For the scale invariant (n = 1) case, Cf^ = [l{l + 1)]"^ 6C|"^. 

^ Example of some inflationary models are: eternal inflation, hybrid inflation, chaotic. Ghost inflation. Tilted Ghost inflation, DBI inflation, brane inflation, 
N-flation, bubble inflation, extended inflation, false vacuum inflation, power law inflation, k-inflation, hyperextended inflation, supersymmetric inflation. 
Quintessential inflation. Natural inflation. Super inflation. Supernatural inflation, D-terni inflation, B -inflation. Thermal inflation, discrete inflation. Assisted 
inflation. Polar cap inflation. Open inflation. Topological inflation. Double inflation. Multiple inflation. Induced-gravity inflation. Warm inflation, stochastic 
inflation. Generalized assisted inflation, self-sustained inflation. Graduated inflation. Local inflation. Singular inflation, SHnky inflation, locked inflation. 
Elastic inflation. Mixed inflation. Phantom inflation. Boundary inflation. Non-commutative inflation, Tachyonic inflation. Tsunami inflation. Lambda inflation, 
Steep inflation. Oscillating inflation. Mutated Hybrid inflation, intermediate inflation, Inhomogeneous inflation. 
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somewhat surprisingly, there exist scenarios where the Universe first contracts and then expands (such as the ekpyrotic/cycHc 
model), which (up to theoretical uncertainties regarding the precise mechanics of the bounce) also reproduce Universes with the 
properties described above. What we would like to do is to find observables that allows us to distinguish between members of the 
inflationary zoo. The exciting fact is that upcoming experiments will have the sensitivity to achieve this goal. Tilt and Running: 
Inflationary models very generically predict a slight deviation from completely flat spectrum. If we write the primordial power 
spectrum as A$(fc) ~ A{ko)(^-^)"^ , then = 1 correspond to flat spectrum and the quantity jris — 1| is called a tilt, which 
characterizes the deviation from scale invariant spectrum. Although the deviations from the scale invariance are predicted to be 
small, the exact amount of deviation depends on the details of the inflationary model. For example in most slow roll models 
|n — 1| is of order l/N^, where A^e ~ 60 is a number of e-folds to the end of inflation. Ghost inflation, however, predicts 
negligible tilt. Hence characterizing the tilt of the scalar spectral index is a useful probe of the early universe. Currently the 
most stringent constraints on tilt come from the WMAP 5-year data, 72s = 0.960l|]:S!i3 ifsill . which already disfavors inflationary 
models with 'blue spectral index' (n^ > 1). The la error on Ug will reduce to Aug — 0.0036 for upcoming Planck satellite and 
to A7is = 0.0016 for futuristic CMBPol like satellite iH. 

Apart from the tilt in the primordial power spectrum, inflationary models also predict 71, to be slightly scale dependent. This 
scale dependence is referred to as 'running' of the spectral index n^, and is defined as drig/d In k. The constraints on the running 
from the WMAP 5-year data are -0.090 < dug/dhik < 0.0019 Etl. The 1 cr error wifl reduce to A{dns/d\nk) = 0.0052 for 
upcoming Planck satellite and to A((ins/dlii A:) = 0.0036 for a fourth-generation satellite such as CMBPol jH]. 

Primordial Gravitational Waves: Inflation also generates tensor perturbations (gravitational waves), which although small 
compared to scalar component, are still detectable, in principle. So far primordial gravitational waves have not been detected. 
There are upper limits on their amplitude; see Ref. ||53I1 for a current observational bounds on the level for primordial gravita- 
tional waves. Detection of these tensor perturbations or primordial gravitational waves is considered a 'smoking gun' for the 
inflationary scenario. In contrast to inflation, ekpyrotic (cyclic) models predict an amount of gravitational waves that is much 
smaller than polarized foreground emission would allow us to see even for an ideal CMB experiment. Primordial scalar pertur- 
bations create only E-modes of the CMB"*, while primordial tensor perturbations generate both parity even E-modes and parity 
odd B-modes polarization [56-58]. The detection of primordial tensor B-modes in the CMB would confirm the existence of 
tensor perturbations in the early universe. This primordial B-mode signal is directly related to the Hubble parameter H during 
inflation, and thus a detection would establish the e nerg y scale at which inflation happened. Various observational efforts are 
underway to detect such B-mode signal of the CMB 1521 . Search for primordial B-modes is challenging. Apart from foreground 
subtraction challenges, and the challenge of reaching the instrumental sensitivity to detect primordial B-modes, there are several 
non-primordial sources such as weak lensing of CMB by the large scale structure I60i,i61il . rotation of the CMB polarization ll62l - 
l65ll . and instrumental systematics that generate B-modes which contaminate the inflationary signal ll66l - l68ll . The amplitude of 
gravitational waves is parametrized as the ratio of the amplitude of tensor and scalar perturbations, r. The limit from WMAP 
5-year data is r < 0.22 (2o-) Isill . 

Isocurvature Modes: Inflationary models with a single scalar field predict primordial perturbations to be adiabatic. Hence 
detection of isocurvature density perturbations is a "smoking gun" for multi-field models of inflation. A large number of infla- 
tionary models with multiple scalar fields predict some amount of isocurvature modes ||69l48l|]. For example, curvaton models 
predict the primordial perturbations to be a mixture of adiabatic and isocurvature perturbations. Isocurvature initial conditions 
specify perturbations in the energy densities of two (or more) species that add up to zero. It does not perturb the spatial curvature 
of comoving slice (i.e TZ is zero, hence the name isocurvature). In general, there can be four types of isocurvature modes, namely: 
baryon isocurvature modes, CDM isocurvature modes, neutrino density isocurvature modes and, neutrino velocity isocurvature 
modes. These perturbations imprint distinct signatures in the CMB temperature and E-polarization anisotropics 18211 . The con- 
tribution of isocurvature modes is model dependent, and different models predict different amounts of it. There exists an upper 
limit on the allowed isocurvature modes using CMB temperature anisotropics Issllsill . a characterization (or detection of any) 
of isocurvature modes has a potential of discriminating between early Universe models. 

Primordial Non-Gaussianity: Canonical inflationary models predict primordials perturbations to be very close to Gaussian isj- 
1, and any non-Gaussianity predicted by the canonical inflation models is very smalljflljll . However models with non- 
linearity ifiolfisl 85], interacting scalar fields |iT2l.^86j|, and deviation from ground state |[87, 88 1 can generate large non-Gaussian 
perturbations. The amplitude of the non-Gaussian contribution to the perturbation is often referred to as /nl even if the nature 
of the non-Gaussianities can be quite different. Different models of inflation predict different amounts of /nl, starting from 
very close to zero for almost Gaussian perturbations, to /nl ~ 100 for large non-Gaussian perturbations. For example, the 
canonical inflation models with slow roll inflation, where only a couple of derivatives of potential are responsible for inflationary 
dynamics, predict /nl ~ 0.05 ifTsll . In models where higher order derivatives of the potential are important, the value of /nl 



To first order in perturbations, primordial scalar perturbations do not generate B-modes of CMB. However at second (and higher) order in perturbations, scalar 
perturbations do produce B-modes IsillsBl . The B-modes generated from higher order perturbations are expected to be smaller than the tensor B-mode levels 
that the upcoming and future experiments (like CMBPol) are sensitive too. 
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folded 

FIG. 3: Shapes of Non-Gaussianity. The shape function F{ki, k2, kz) forms a triangle in Fourier space. The triangles are parametrized by 
rescaled Fourier modes, 12 ~ k2/k\ and 13 = k^/ki. Figure from Ref. (43I1 

varies from /nl ^0.1 where higher order derivatives are suppressed by a low UV cutoff ll89ll to /nl ^ 100 based on Dirac- 
Born-Infeld effective action. Ghost inflation, where during inflation, the background has a constant rate of change as opposed 
to the constant background in conventional inflation, is also capable of givin g /n t, ^ 100 190,1. The additional field models 
generating inhomogeneities in non-thermal species ll9lll can generate /nl ^ 5 ll92ll : while curvaton models, where isocurvature 
perturbations in second field during the inflation generate adiabatic perturbations after the inflation, can have /nl '--^ 10 ll93ll . 

In the following we will see that non-Gaussianity, far from being merely a test of standard inflation, may reveal detailed 
information about the state and physics of the very early Universe, if it is present at the level suggested by the theoretical 
arguments above. 

in. PRIMORDIAL NON-GAUSSIANITY 

Large primordial non-Gaussianity can be generated if any of the following condition is violated ll94l] 

• Single Field. Only one scalar field is responsible for driving the inflation and the quantum fluctuations in the same field is 
responsible for generating the seed classical perturbations. 

• Canonical Kinetic Energy. The kinetic energy of the field is such that the perturbations travel at the speed of light. 

• Slow Roll. During inflation phase the field evolves much slowly than the Hubble time during inflation. 

• Initial Vacuum State. The quantum field was in the Bunch-Davies vacuum state before the quantum fluctuation were 
generated. 

To characterize the non-Gaussianity one has to consider the higher order moments beyond two-point function, which contains 
all the information for Gaussian perturbations. The 3-point function which is zero for Gaussian perturbations contains the infor- 
mation about non-Gaussianity. The 3-point correlation function of Bardeen's curvature perturbations, can be simplified 
using the translational symmetry to give 

($(ki)$(k2)$(k3)) = (27r)353(ki + k2 + k3)/NL • F{ki,k2, h) . (17) 

where F{ki, ^2,^3) tells the shape of the bispectrum in momentum space while the amplitude of non-Gaussianity is captured 
dimensionless non-linearity parameter /nl- The shape function F{ki,k2, fcs) correlates fluctuations with three wave-vectors 
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FIG. 4: Plot of the function -F(l, X2, X3)x2X3 for the Slow-Roll inflation as given by Eq. (21) (left panel) and the local distribution as given by 
Eq. (19) (right panel). The figures are normalized to have value 1 for equilateral configurations x2 = x3 = 1 and set to zero outside the region 
1 — X2 < x-^ < X2- Here xZ = fc3/fcl, x2 = fc2/fcl and (. = r) — 1/30. The figures are taken fromBabich et al. 2003 (9^ . 



and form a triangle in Fourier space. Depending on the physical mechanism responsible for the bispectrum, the shape of the 
3-point function, F{ki, k2,k^) can be broadly classified into three classes ll95ll . The local, "squeezed," non-Gaussianity where 
F{ki, k2, fca) is large for the configurations in which fci «C fc 2 ~ k^. Most of the studied inflationary and Ekpyrotic models 
produce non- Gaussianity of local shape (eg. j?! IgllgeHTTl). Second, the non -local, "e quilateral," non-Gaussianity where 
F{ki , fc2, fca) is large for the configuration when fci « ^2 ~ fcs. Finally the folded ll l4ll isll shape where F{ki , k2, fcs) is large 
for the configurations in which fci « 2k2 ~ ik^. Figure |3] shows these three shapes. 

Non-Gaussianity of Local Type: The local form of non-Gaussianity may be parametrized in real space as^ ifislll 16lll ITll : 



C(r) = CL(r) + -/NL (C!(r)-(C!(r))) 



(18) 



where Cl(^) is the linear Gaussian part of the perturbations, and /nl characterizes the amplitude of primordial non-Gaussianity. 
Different inflationary models predict different amounts of /nl, starting from 0(1) to /nl ^ 100, beyond which values have 
been excluded by the Cosmic Microwave Background (CMB) bispectrum of WMAP temperature data. The bispectrum in this 
model can be written as 



Fiocaiih, k2,h) = 2A|/nl 



1 



1 



1 



,3-(n,-l),3-(n,-l) 
'^1 '^2 



-l),3-(n, 
«,3 



, 3-(n, 



-l)iL3- 



(19) 



where A$ is the amplitude of the primordial power spectrum. 

The local form arises from a non-lin ear relatio n between inflaton and curvature perturbations curvaton mod- 

els or the New Ekpyrotic models Ill8llll9ll . Models with fluctuations in the reheating efficiency [9, 10] and multi-field 
inflationary models [17] also generate non-Gaussianity of local type. 

Being local in real space, non-Gaussianity of local type describes correlations among Fourier modes of very different k. In 
the limit in which one of the modes becomes of very long wavelength |l2Cl|| . k^ -> 0, (i.e. the other two fc's become equal 
and opposite), freezes out much before fci and ^2 and behaves as a background for their evolution. In this limit Fiocai is 
proportional to the power spectrum of the short and long wavelength modes 



-Piocal OC 



1 1 



As an example, for canonical single field slow-roll inflationary models, the three point function is given by 11151 



Fsiow-miiiki,k2,k3) = i(3/5)A|j^ 



(3e 



kf 



8e 



'fci 



(20) 



(21) 



5 Orequivalently $(r) = '5>L{r) + /nl (*i(r) - (^lir))) 
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F(X2, X3) 



FIG. 5: Plot of the function F{1, X2, X3)x2x'^ for the inflation with higher derivatives as given by Eq. (23) (left panel) and the ghost inflation 
as given by Eq. (24) (right panel). The figures are normalized to have value 1 for equilateral configurations x2 = x3 = 1 and set to zero outside 
the region 1 — X2 < x^ < X2- Here x3 = k3/kl and x2 = fc2/fcl. The figures are taken from Babich et al. 2003 (9^ 



where e and 77 are the usual slow-roll parameters and are assumed to be much smaller than unity. Taking the limit — > 
gives the local form as in Eq. ( |20l ). To show this point, Figure. |4] compares the non-Gaussianity shape for local type and for 
slow-roll model. Although in this limit, slow-roll models do predict no-Gaussianity of local type but as evident from Eq. ( I2TI ). 
the bispectrum of inflaton perturbations yields a non-trivial scale dependence of /nl lfl2l[l5tl . However in the slow roll limit 
?7, e << 1 and hence the amplitude is too small to detect. 

Non-Gaussianity of equilateral type: While vast number of inflationary models predict non-Gaussianity of local type, this 
model, for instance, fails completely when non-Gaussianity is localized i n a s pecific range in k space, the case that is predicted 
from inflation models with higher derivative kinetic terms ll90llll5lll2l[4l24ll . In these models the correlation is among modes 
with comparable wavelengths which go out of the horizon nearly at the same time. The shape function for the equilateral shape 
can be written as ll25ll 



1 



,3-(n,-l), 3-(«,-l) 
"-1 "'2 



1 



l-(n,-l)/3,^2-2(ri,-l)/3^3-(ri,-l) 



(2 perm.) 



(5 perm. 



(fclfc2fc3)2-2(».-l)/3 



(22) 



The models of this kind have large F{ki, ^2, k^) for the configurations wh ere fc i « ^2 ~ ^3. The equilateral form arises 
from non-canonical kinetic terms such as the Dirac-Born-Infeld (DBI) acti on 1 122 1. the ghost condensation llooll . or any other 
single-field models in which the scalar field acquires a low speed of sound ll IsT 124ll . 

As an example, models with higher derivative operators in the usual inflation scenario and a model of inflation based on the 
Dirac-Born-Infeld (DBI) action produce a bispectrum of the form 



F,,{kuk.,k3) = ^Al^^^ 



1 



ki + k2 + k^ 



(23) 



The above model uses g^(V(/))2(V(/))2 as a leading order operator. DBI inflation, which can produce large non-Gaussianity, 
/nl ^ 100, also has F{ki, k2, k^) of a similar form. 

Ghost inflation, where an inflationary de Sitter phase is obtained with a ghost condensate, produces a bispectrum of the 
following form ll90ll 



Fghost{ki,k2,k3) = -(3/5)' 
1 

TPF 



' r(l/4)3 Af2" 
di]i]-^F*{i])F* (^,^]F'* 



(24) 



k^' 



V ] ^3(^1 ■ ^2) + symm., 
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where a and /3 are free parameters of order unity, and 

F{x) = y|(^x)^/^i/^;i(xV2). (25) 

Ghost inflation also produces large non-Gaussianity, /nl ~ 100. Figure shows the shape of non-Gaussianity of equilateral 
type by showing F{ki, ^2,^3) for ghost inflation and for a model with a higher derivative term. 

Folded Shape: So far the 3-point functions were calculated assuming the regular Bunch-Davis vacuum state, giving rise to 
either local or equilateral type non-Gaussianity. However if the bispectrum is calculated by dropping the assumption of Bunch- 
Da yis initial state give rise to bispectrum shape which peaks for the folded shape, ki w 2fc2 ~ 2fc3, with shape function given 

asiinillillill 

F„™.BD(fci,fc2,fc3) = (3/5)' • _i-^7^e(/3fe^) (cos(fc,^o) - l) (26) 

where /Jfe^ are the Bogoliubov coefficients which encode information about the initial conditions, 770 is the initial conformal time 
and kj = ~ 2fci. 



IV. THE COSMIC MICROWAVE BACKGROUND BISPECTRUM 



Since the discovery of CMB by Penzias and Wilson in 1965 Hi 2611 and the first detection of CMB temperature anisotropies on 
large scales by the COBE DMR lll27ll . the space satelUte WMAP and over a dozens of balloon and ground based experiments 
have characterized the CMB temperature anisotropies to a high accuracy and over a wide range of angular scales. The space 
satellite Planck which launched in 2009 will soon characterize the temperature anisotropies to even higher accuracy up to angular 
scales of imax ~ 2500. The CMB power spectrum is obtained by reducing all the information of Npi^^ (^ 10^ for WMAP and 
^ 10^ for Planck). Such reduction is justified to obtain a fiducial model, given the non-Gaussianities are expected to be small. 
With high quality data on the way, the field of non-Gaussianity is taking off. CMB bispectrum contains information which is not 
present in the power-spectrum and as we say in the previous section, is a unique probe of the early universe. 

The harmonic coefficients of the CMB anisotropy aim = j '^^nAT(n)y^* ^ can be related to the primordial fluctuation 
$ as 

/d^k 
^^k)g^,{k)Y*Jk) + n,m, (27) 

where ^(k) is the primordial curvature perturbations, for a comoving wavevector k, g^{r) is the radiation transfer function 
where the index p refers to either temperature (T) or /s-polarization (E) of the CMB. A beam function and the harmonic 
coefficient of noise nem are instrumental effects. Eq. (l27T i is written for a flat background, but can easily be generalized. 

Any non-Gaussianity present in the primordial perturbations $(k) gets transferred to the observed CMB via Eq. (l27l ). The 
most common way to look for non-Gaussianity in the CMB is to study the bispectrum, the three-point function of temperature 
and polarization anisotropies in harmonic space. The CMB angular bispectrum is defined as 

fjPqr _ I p q r \ /^ox 

-"£i£2fe,mim2m3 — \"f imi "f2m2 "fsms ' ' ^ °' 

and the angular-averaged bispectrum is 

^pqr _ ^ 1^1 ^2 \ T^pqr 



jypqi _ X -J jypqi 



where the matrix is the Wigner 3J symbol imposing selection rules which makes bispectrum zero unless 

(i) ^1 + £2 + h ^integer 

(ii) mi + mi + 7713 = 

(iii) \l, -F^<ik<h + £j for i,j,k = l, 2, 3. 
Using Eq.dZTli the bispectrum can be written as 

xgl{ki)gl{k2)glik,) (<i>(ki)$(k2)<I>(k3)) , (30) 
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where ($(ki)$(k2)$(k3)) is the primordial curvature three-point function as defined in Eq. ( fTTI ). 

To forecast constraints on non-Gaussianit y usi ng CMB data, we will perform a Fisher matrix analysis. The Fisher matrix for 
the parameters pa can be written as 



J'ab 



E E ■ 

{ijk, pqr} ii<i2^^3 



dpa hjk,pqr Qp^ 



(31) 



The indices a and b run over all the parameters bispectrum depends on, we will assume all the cosmological parameters except 
/nl to be known. Indices ijk and pqr run over all the eight possible ordered combinations of temperature and polarization 
given by TTT, TTE, TET, ETT, TEE, ETE, EET and EEE; the combinatorial factor A^.^^^a equals 1 when all ts are 
different, 6 when £1=^2= ^3, and 2 otherwise. The covariance matrix Cov is obtained in terms of Cf^, Cf^, and Cj^ 
(see ll20l[34ll ) by applying Wick's theorem. 

For non-Gaussianity of the local type, for which the functional form F{ki, k2, k^) is given by Eq. (fT9] l. we have 



8B "-^^ 



dfNL V 47r 
where the functions a and /3 are given by 



(2^i + l)(2£2 + l)(2£3 + l) f h h h 





2 ; r 
'0 



-al^l3ll3l+ 2 perm. 



(32) 



alir) 



- jdk e glik) Mkr) , 

TT 



dk k-^ gl(k) jiikr)A^ /fc"'"^ , 



(33) 
(34) 



In the expression above we use the dimensionless power spectrum amplitude A$, which is defined by P$(fc) = A^k 1)^ 
where Ug is the tilt of the primordial power spectr um. O ne can compute the transfer functions gj{k) and gf{k) using publicly 
available codes such as CMBfast IH and CAMB lfl28ll 

In a similar way, from Eq. (|22] |. one can derive the following expressions for the bispectrum derivatives in the equilateral case. 



dB 



ijk 



dfi 



NL 



(2£i + l)(2£2 + l)(2£3 + l) f ei £2 £3 





4tt 



X 6 jr^dr 
where the functions 6, and 7 are given by 



-<PiPt, + 2 perm. + PlllS^ + ^ PC^^^^- " 2<5^/^/, 



(35) 



^ - dkkg\{k)3i{kr)A'J^k^-^-^^/\ 
2 



TT 



dkgj{k) ji{kr)A 



2/3 , 2(„,-l)/3 



(36) 
(37) 



Recently a new bispectrum template shape, an orthogonal shape has been introduced 1112911 which characterizes the size of the 
signal if^l^") which peaks both for equilateral and flat-triangle configurations. The shape of non-Gaussianities associated with 
fwL^° is orthogonal to the one associated to /^'^'' . The bispectrum for orthogonal shape can be written as lll2' 



dB 



ijk 



ortho 
NL 



(2£i + l)(2^2 + l)(2£3 + l) ( tl ^2 ^3 





47r 



X 18 



'O^iAA. + 2 perm. + Pl.ll^Sl^ + 5 perm. - -5\^5ljl 



(38) 



Estimator 



An unbiased bispectrum-based minimum variance estimator for the nonlinearity parameter in the limit of full sky and homo- 
geneous noise can be written as UtI I21II25I1 



/nl = 



N ^ 



Bi 



H TO2 m3 J Ci^Ce^d 



^iimi 0'i2m2 ^^3?Tl3 



(39) 
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where Bi-^i^i^ is angle averaged theoretical CMB bispectrum for the model in consideration. The normalization N can be 
calculated to require the estimator to be unbiased, (/nl) = /nl- If the bispectrum Bp^i^i^ is calculated for /nl = 1 then the 
normaUzation takes the following form 

^-Et^T^ (40) 
The estimator for non-Gaussianity, Eq. ( [39] l. can be simplified using Eq. jZTl l to yield 

oo 

•^NL " ' E / d^fiYi^,nAn)yhinAn)Yi.,„r.,{h) j r^dr ji,{kir)ji^{k2r)ji.,{k3r) C-^C-^C-^ 



liTUi 

2kldki 2kldk2 2k^dk 



■F{ki,k2,k3)Af_^{ki)Af^{k2)Af^{k3) ai^rinai^m^ai^rna , (41) 



where F{ki, k2, fcs) is a shape of 3-point function as defined in Eq. iT% . Given the shape F{ki, k2, fcs), one is interested in, 
it is conceptually straightforward to constrain the non-linearity parameter from the CMB data. Unfortunately the computation 

time for the estimate scales as Np(^, which is computationally challenging as even for the WMAP data the number of pixels is 

fi 3/2 

of order Npi^ ^ 10 . The scaling can be understood by noting that the each spherical harmonic transform scales as N^^^ and 
the estimator requires £^(a Npix) number of spherical harmonic transforms. 
The computational cost decreases if the shape can be factorized as 

F{ki,k2, fcs) = h{h)h[k2)h{k3) , (42) 

with which the estimator simplifies to 

./nl = ^-Jd^h J r^dr he/ ^^Jk{kr)Mk)Al{k)Cr^'ai^mri^rr.M (43) 





3/2 7 

and computational cost now scales as A^^/^ . For Planck {Npi^ ^ 5 x 10 ) this translates into a speed-up by factors of millions, 
reducing the required computing time from thousands of years to just hours and thus making /nl estimation feasible for future 
surveys. The speed of the estimator now allows sufficient number of Monte Carlo simulations to characterize its statistical 
properties in the presence of real world issues such as instrumental effects, partial sky coverage, and foreground contamination. 
Using the Monte Carlo simulations it has been shown that estimator is indeed optimal, where optimality is defined by saturation 
of the Cramer Rao bound, if noise is homogeneous. Note that even for the non-factorizable shapes, by usi ng th e flat sky 
approximation and interpolating between the modes, one can estim ating /nl in a computationally efficient way ll30ll . 

The extension of the estimator of /nl from the temperature data lUTIl to include both the temperature and polarization data of 
the CMB is discussed in Yadav et al. Summarizing briefly, we construct a cubic statistic as a combination of 

(appropriately filtered) temperature and polarization maps which is specifically sensitive to the primordial perturbations. This 
is done by reconstructing a map of primordial perturbations, and using that to define a fast estimator We also show that this 
fast estimator is equivalent to the optimal estimator by demonstrating that the inverse of the covariance matrix for the optimal 
estimator lt34ll is the same as the product of inverses we get in the fast estimator The estimator still takes only N^^^ operations 

5 /2 

in comparison to the full bispectrum calculation which takes N^^l^ operations. 

S 

For a given shape the estimator for non-linearity parameter can be written as /nl = „° , where for the equilateral, local 

shape 

and orthogonal shapes, the Sshape can be written as 

SequUaterai = / dr [ d^ h[B {h, r)B (h, r)Aih, r) + ^Din, r)^ - 2B {fi, r)C (h, r)Din, r)] (44) 

J sky J J 

Siocal - 7^ fr^dr [ d^nB{n,r)B{h,r)A{h,r) (45) 

Jsky J J 

Sorthogonai = J j d^ h[B {h, r)B {h, r)Aih, r) + ^D{n, - 2B {h, r)C (h, r)Din, r)] (46) 

with 

ip Im ip Im 
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ip hn ip hn 

and fsky is a fraction of sky. Index i and p can either be T or E. 

^= E E A^^"™(^"')^^(^"')i'(^"')t^0^r' (49) 

ijkpqr 2<li<l2<l3 ^i^a^a 

Indices 1, 7, k,p, q and r can either be T or iJ. Here, A^j^^^g is 1 when ii £2 ^ £3, 6 when ^1 = £2 = ^3, and 2 otherwise, 
^fXf"™ '-^^ theoretical bispectrum for /nl = 1 llzoll . 

It has been shown that the above estimators defined in Eq. ( |46] | are minimum variance amongst bispectrum-based estimators 
for full sky coverage and homogeneous noise ll20ll . To be able to deal with the realistic data, the estimator has to be able to deal 
with the inhomogeneous noise and foreground masks. The estimator can be generalized to deal with partial sky coverage as well 
as inhomogeneous noise by adding a linear term to Sprim- Sprim Sprim + S^priin^ ■ For the temperature only case, this has 
been done in ll25ll . Following the same argument, we find that the linear term for the combined analysis of CMB temperature 
and polarization data is given by 

Spr^'' = "T^ / ^''^^ / d^f" {25("' r))MC + A{h, r) {B%^{n, r))Mc} , (50) 

J sky J J 

where Asim{n, r) and Bsim{fi, r) are the A and B maps generated from Monte Carlo simulations that contain signal and noise, 
and (..) denotes the average over the Monte Carlo simulations. 
The generalized estimator is given by 

C" I olinear 

~ '-'prim -r J prim 

./NL = j;^ • Pl) 

Note that {Sp^^"^) mc — ~{Sprim) mc, ™d this relation also holds for the equilateral shape. Therefore, it is straightforward to 
find the generalized estimator for the equilateral shape: first, find the cubic estimator of the equilateral shape, SequU., and take 
the Monte Carlo average, (SequU.) mc- Let us suppose that SequU. contains terms in the form of ABC, where A, B, and C are 
some filtered maps. Use the Wick's theorem to re-write the average of a cubic product as {ABC)mc = {A) mc{BC) mc + 
{B)mc{AC)mc + (C) MC (AB) MC ■ Finally, remove the MC average from single maps, and replace maps in the product with 
the simulated maps {A)mc{BC)mc + {B)mc{AC)mc + {C)mc{AB)mc ^ A{BsimCsim)MC + B{As.imCsim)MC + 
C {AsimB Sim) MC ■ This operation gives the correct expression for the linear term, both for the local form and the equilateral 
form. 

The main contribution to the linear term comes from the inhomogeneous noise and sky cut. For the temperature only case, 
most of the contribution to the linear term comes from the inhomogeneous noise, and the partial sky coverage does not contribute 
much to the linear term. This is because the sky-cut induces a monopole contribution outside the mask. In the analysis one 
subtracts the monopole from outside the mask before measuring Sprim, which makes the linear contribution from the mask 
small 12511 . For a combined analysis of the temperature and polarization maps, however, the linear term does get a significant 
contribution from a partial sky coverage. Subtraction of the monopole outside of the mask is of no help for polarization, as the 

monopole does not exist in the polarization maps by definition. (The lowest relevant multipole for polarization is I = 2.) 

3/2 

The estimator is still computationally efficient, taking only (times the r sampling, which is of order 100) operations in 
comparison to the full bispectrum calculation which takes N^^^ operations. Here Npi^ refers to the total number of pixels. For 
Planck, Npix ^ 5 x 10^, and so the full bispectrum analysis is not feasible while our analysis is. 



V. CONSTRAINTS FROM THE CMB BISPECTRUM 



A. Current Status 



Currently the the Wilkinson Microwave Anisotropy Probe (WMAP) satellite provides the "best" (largest number of signal 
dominated modes) CMB data for non-Gaussianity analysis. Over the course of WMAP operation the field of non-Gaussianity 
has made vast progress both in terms of theoretical predictions of non-Gaussianities from inflation and improvement in the 
bispectrum based estimators. At the time of WMAP's first data release in 2003 the estimator was sub-optimal in the presence of 
partial sky coverage and/or inhomogeneous noise. With the sub-optimal estimator one could not use the entirety of WMAP data 
and only the data up to Imax 350 were used to obtain the constraint = 38 ± 96(2cr) ||23f|. These limits were around 30 

times better than the previous constraints of |/nl | < 1.5 x 10'^ from the Cosmic Background Explorer (COBE) satellite [? ]. 
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By the time of second WMAP release in 2007 the estimator was generalized by adding a linear to the KSW estimator which 
allows to deal with partial sky coverage and inhomogeneous noise. The idea of adding a linear term to reduce excess variance 
due to noise inhomogeneity was introduced in ll25ll . Apphed to a combination of the Q, V and W channels of the WMAP 3 -year 
data up to l,nax ^ 400 this estimator had yielded the tightest constraint at the time on /nl as: —36 < /nl < 100(2(t) i26l| . This 
estimator was further generalized to utilize both the temperature and E-polarization information in ll2lll . where it was pointed 
out that the linear term had been incorrectly implemented in Eq. 30 of 12511 ■ Using Monte-Carlo simulations it has been shown 
that this corrected estimator is nearly optimal and enables analysis of the entire WMAP data without suffering from a blow-up 
in the variance at high multipoles*. The first analysis using this estimator shows an evidence of non-Gaussianity of local type 
at around 2.8cr in the WMAP 3-year data. Independent analysis shows the evidence of non-Gaussianity at lower significance, 
around 2.5(7 (see TablelVlt. 

By the time of the third WMAP data release (with 5-year obsevational data) in 2008 the /nl estimation technique was 
improve d fur ther by implementing the covariance matrix including inhomogeneoous noise to make the estimator completely 
optimal II131II . Using the optimal estimator and using the entire ty of WMAP 3-year data there is an evidence for non-Gaussianity 
of local type at around 2.5a level /nl ~ 58 ± 23(1(t) |I131[1 . However with WMAP 5-year data the significance goes down 
from ~ 2.5ct to ~ I.Sct |I131I1 . Table IVIIcompares the constraints obtained by different groups using WMAP 3-year and WMAP 
5-year data. Fig|6]shows this comparison in more detail, showing the constraints also as a function of maximum multipole Imax 
used in the analysis. Few comments are in place: 1) constraints on /nl from WMAP 3-year data as a function of Imax show 
a trend where the mean value rises at around Imax = 450, below which data is consistent with Gaussianity and above which 
there is deviation from Gaussianity at above 2fT. The result becomes roughly independent of Imax > 550 with evidence for 
non-Gaussianity at around 2.5(T level, 2) independent analysis and using different estimators (optimal and near-optimal with 
linear term) sees this deviation from non-Gaussianity at around 2.5<t in WMAP 3-year data, 3) significance of non-Gaussianity 
goes down to around 2a with WMAP 5-year data. The drop in the mean value between WMAP 3-year and 5-year data can be 
attributed to statistical shift. 

The best constraints on the equilateral and orthogonal shape of non-Gaussianity using the WMAP 5-year data are f^^'^ = 
155 ± 140(lcr) and /°'l°^- = 149 ± llO(lcr) respectively lfl2i . 

As we were completing this article, the WMAP 7-year data was r eleased, with constraints f^l^^ = 32 ± 21(1ct), /^7'- = 
26 ± 140(lcr) and /°'l°^ -202 ± 104(lcr) fl32ll . 



B. Future Prospects 

Now we discuss the future prospects of using the bispectrum estimators for constraining the non-linearity parameter /nl 
for local and equilateral shapes. We compute the Fisher bounds for three experimental setups, (1) cosmic variance limited 
experiment with perfect beam (ideal experiment hereafter), (2) Planck satellite with and noise sensitivity Ap = 56/iK-arcmin 
and beam FWHM a = T (3) a futuristic CMBPol like satellite experiment with noise sensitivity Ap = 1.4/iK-arcmin and 
beam FWHM cr = 4' (CMBPol hereafter). Beside /nl we fix all the other cosmological parameters to a standard fiducial 
model with a flat ACDM cosmology, with parameters described by the best fit to WMAP 5-year results jsi]], given by ilb = 
0.045, flc = 0.23, Hq = 70.5, Us = 0.96, rit = 0.0, and t = 0.08. We calculate the theoretical CMB ti-ansfer functions and 
power spectrum from publicly available code CMB FAST ifsoll . We also neglect any non-Gaussianity which can be generated 
during recombination or there after We discuss the importance and effect of these non-primordial non-Gaussianities in next 
section. 

The scaUng of signal-to-noise as a maximum multipole Imax for the local Il33lll34ll and equilateral model lll35ll are 

{S/N)iocal OC hvtmax , {S/N)equil OC \/ tmax ■ (52) 

In principle one could go to arbitrary high Imax but in reality secondary signals will certainly overwhelm primary signal beyond 
Imax > 3000, we restrict to the analysis to Imax — 3000. In Figure |7] we show the la Fisher bound as function of maximum 
multipole imax, for local and equilateral type of non-Gaussianity. For both local and equilateral case we show the Fisher bound 
for the analysis using only the CMB temperature information (TTT), only the CMB polarization information (FEE), and the 
combined temperature and polarization analysis. Note that by having both the temperature and £'-polarization information one 
can improve the sensitivity by combining the information. Apart from combining the T and E signal, one can also do cross- 
checks and diagnostics by independently analysing the data. Temperature and polarization will have different foregrounds and 
instrumental systematics. 



We will refer to the estimator in Ref. as a near-optimal-vl, while the corrected estimator of Ref. (2^ as near-optimal. 
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TABLE I: Hint of local non-Gaussianity at 2g level 



data (mask, estimator) 


/nT"' ± la error 


deviation from Gaussianity 




WMAP 3-year (KpO, near-optimal) 


87 ±31 


2.8(T 


Yadav and Wandelt [28J 


WMAP 3-year (KQ75, optimal) 


58 ±23 


2.5cr 


Smith et al. [1311 


WMAP 3-year (KpO, near-optimal) 


69 ±30 


2.3cr 


Smith et al. [1311 


WMAP 5-year (KQ75, near-optimal) 


51 ±30 


1.7cr 


Komatsu et al. [511 


WMAP 5-year (KQ75, optimal) 


38 ±21 


l.Scr 


Smith et al. [1311 



The difference between the results by Ref. and Il3lll for WMAP 3-year data using the KpO mask can be a result of different choices of weighting in the 
near-optimal estimator. The optimal estimator has a unique weighting scheme. 



A CMBPol like experiment will be able to achieve the sensitivity of A/^"]^"' ~ 2(1(t) for non-Gaussianity of local type 
and A/X - 13(lcr) for non-Gaussianity of equilateral type. For the local type of non-Gaussianity this amounts to an 
improvement of about a factor of 2 over the Planck satellite and about a factor of 12 over current best constraints. These 
estimates assume that foreground cleaning can be done perfectly, i.e. the effect of residual foregrounds has been neglected. 
Also the contribution from unresolved point sources and secondary anisotropics such as ISW-lensing and SZ-lensing has been 
ignored. 

Running non-Gaussianity: The primordial non-Gaussian parameter /nl has been shown to be scale-dependent in several 
models of inflation with a variable speed of sound, such as Dirac-Born-Infeld (DBI) models. Starting from a simple ansatz for a 
scale-dependent amplitude of flie primordial curvature bispectrum for primordial non-Gaussianity, 



/nl ^ ./nl t- (53) 




where K = (fcifc2fc3)^^'^ and kp is a pivot point. The primordial bispectrum is therefore determined in terms of two parameters: 
the amplitude /nl and the new parameter tinc quantifying its running. One can generalize the Fisher matrix analysis of the 
bispectra of the temperature and polarization of the CMB radiation and derive the expected constraints on the parameter n 
that quantifies the running of /NL(fc) for current and future CMB missions such as WMAP, Planck and CMBPol. We will 
consider some non-zero /nl as our fiducial value for the Fisher matrix evaluation. Clearly, in order to be able to constrain a 
scale-dependence of /nl, its amplitude must be large enough to produce a detection. If /nl is too small to be detected (/nl < 2 
is a lowest theoretical limit even for the ideal experiment), we will obviously not be able to measure any of its features, either. 
In the following we will then always consider a fiducial value of /nl large enough to enable a detection. Figure [8] shows the 
I — (7 joint constraints on /nl and ung- In the event of a significant detection of the non- Gaussian component, corresponding 
to /nl = 50 for the local model and /nl = 100 for the equilateral model of non-Gaussianity, is able to determine ting with 
a 1 — (7 uncertainty of An^G — 0.1 and Ati^vg — 0.3, respectively, for the Planck mission and a factor of two better for 
CMBPol. In addition to CMB one can include the information of the galaxy power spectrum, galax y bis pectrum, and cluster 
number counts as a probe of non- Gaussianity on small scales to further constrain the two parameters |l 36|l . 



C. Contaminations 



A detection of non-Gaussianity has profound implications on our understanding of the early Universe. Hence it is impor- 
tant to know and quantify all the possible sources of non-Gaussianities in the CMB. Here we highlight some sources of non- 
Gaussianities due to second-order anisotropics after last scattering surface and during recombination. The fact that Gaussian 
initial conditions imply Gaussianity of the CMB is only true at linear order. We will also discuss the effects of instrumental 
effects and uncertainties in the cosmological parameters on the bispectrum estimate. 



Secondary non-Gaussianities 

Current analysis of the CMB data ignore the contributions from the secondary non-Gaussianities. For WMAP resolution 
it may not be a bad approximation. Studies o f the dominant secondary anisotropics conclude that they are negligible for the 
analysis of the WMAP data for l„,ax < 800 ITl7i[l37ll . However on smaller angular scales several effects starts to kick in, 
for example, 1) the bispectrum contribution due to unresolved point source like thermal Sunyaev-Zeldovich clusters or standard 
radio sources, 2) three way correlations between primary CMB, lensed CMB and secondary anisotropics. We will refer to the 
bispectrum generated due to this three way correlation as B^'^condar y-K wheresome secondaries are, the integrated Sachs-Wolfe 
(ISW) B^'^^'-", Sunyaev-Zeldovich signal and Rees-Sciama H [YTtI [T37ijl40ll . 
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FIG. 6: Top panel: Constraints on local /nl using WMAP 3-year data as a function of maximum multipole Imax used in the analysis. The 
red circles are the results obta ined using near optimal estimator by Yadav & Wandelt f2§l\ . The green triangles are using the the near-optimal 
estimator by Smith et al. Il3lh . The blue square results are obtained using either the optimal estimator by Smith et al. j 13 ill . For all the three 
analysis KpO mask was used. Bottom panel: Comparison between 5-year resul ts (op timal estimator, raw maps) reported in Komatsu et al. (Slll 
and results obtained using the optimal or suboptimal estimator by Smith et al. Il3lll . 



For Future experiments such as Planck and CMBPOl the joint estimation of primordial and secondary bispectrum will be 
required. The observed bispectrum in general would take the following form 

^iiLis = /nl^£i72™3 + bp^B^^i^is + ^szB^fi'i'^ + AjswBg^^g~'^ + ... (54) 

The amplitude of bispectrum due to primary-lensing-secondary cross-correlation is proportional to the product of primary CMB 
power-spectrum and power spectrum of cross-correlation between secondary and lensing signal. 

The reduced bispectrum from the residual point sources (assuming Poisson distributed) is constant i.e. b^^i^i^ = constant. The 
value of the constant will depend on the flux limit at which the point source can be detected and on assumed flux and frequency 
distribution of the sources. 
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FIG. 7: Fisher predictions for minimum detectable /nl (at 1 cr) as a function of maximum multipole Imax- Upper panels are for the local 
model while lower panels are for the equilateral model. Left panels shows an ideal experiment, middle panels are for CMBPol like experiment 
with noise sensitivity Ap = 1.4/^K-arcmin and beam FWHM a — 4' and right panels are for Planck like satellite with and noise sensitivity 
Ap — 40^K-arcmin and beam FWHM cr = 5'. In all the panels, the solid lines represent temperature and polarization combined analysis; 
dashed lines represent temperature only analysis; dot-dashed lines represent polarization only analysis. 



Depending on the shape of primordial bispectrum in consideration, some secondary bispectra are more dangerous than others. 
For example, ISW-lensing B^^^~'^ peaks at the "local" configurations, hence is more dangerous for local primordial shape 
than the equilateral primordial shape. For example for the Planck satellite if the secondary bispectrum is not i ncorp orated 
in the analysis, the ISW-lensing contribution will bias the estimate for the local /nl by around A/^"™' « 10 114111 . The 
bispectrum contribution from primary-lensing-Rees-Sciama signal also peaks at squeezed limit and contribute to effective local 
j-iocai^ 10 I142II . For Planck sensitivity the point source will contamination the l ocal non-Gaussianity by around A/^°l°' ^ 
1 I143II . A recent analysis of the full second order Boltzmann equation for photons ll44ll claims that second order effects add a 
contamination A/nl ^ 5. 

The generalization of the Fisher matrix given by Eq. ( |3TI ) to include multiple bispectrum contribution is 

•^afa - 2. 2. ^ g , (55) 

{gk, pqr}li<l2<ia 

where the additional X and Y indices denote a component such as primordial, point-sources, ISW-lensing etc. For fixed cosmo- 
logical parameters, the signal-to-noise {S/N)i for the component i is 

^ (56) 




Non-Gaussianities from recombination 



Non-Gaussianities can be generated during recombination. One requires to solve secon d order Bo ltzmann a nd Ei nstein equa- 
tions to evaluate the effect. The second order effect on CMB is an active field of study lf54lll45l4l63ll . see Ref. Il64ll for a recent 
review. The non-Gaussianities produced during recombination comprise of various effects, for example see Ref. ll54ll . 

The dominant bispectrum due to perturbative recombination comes from perturbations in the electron density. The amplitude 
of perturbations of the free electron density 5e is around a factor of 5 larger than the baryon density perturbations l\6y\. The 
bispectrum gene rated due to 5e peaks around the "local" configuration with coiTesponding effective non-linearity amplitude 
/nl - few ilEinEi. 



18 



0.4 



0.2 



I 0.0 



-0.2 



-0.4 



Planck 
/„„=2500 



0.5 



I 0.0 



-0.5 




30 40 



50 

/nl 



60 70 



50 100 150 200 

/ml 



0.2 



0.1 



I 0.0 



-0.1 



-0.2 



CMBPol 

Amx=3000 




: I 

■ 1 



Planck 








500/ 












a 


3)f 

✓ 
















FIG. 8: l-a constraints on /nl and ting for local (upper panels) and equilateral model (lower panels) assuming kp — 0.04 Mpc^ and fiducial 
values /nl = 50, hng ~ for the local case and /nl = 100, ting = for the equilateral case. Dashed lines correspond to the limits from 
the temperature information alone, dotted lines to polarization (EEE), while the continuous lines correspond to all bispectrum combinations. 
We consider Planck (left panels) and CMBPOl like (right panels) CMB experiment. 



The b ispec trum contribution due to second order terms which are the products of the first order perturbations is calculated 
in Ref. il69r . The bispectrum contribution from these terms which also peak for the squeezed triangles is small and can be 
neglected in the analysis. For example the signal-to-noise is about 0.4 at Imax = 2000 for a full-sky, cosmic variance limited 
experiment. 

Another contribution to bispectrum which peaks for the equilateral configurations comes from the non-linear evolution of the 
second order gravitational potential. Because of this effect the minimum detectable non-Gaussianity parameter /^l"*' changes 
by A/^l"*'' = 0(10) for Planck like experiment lll35ll . The bispectrum peaks for the equilateral shape because the the growth 
of potential happens on scales smaller than the horizon size. 

On large scales, in the absence of primordial n on-G aussianities and assuming matter domination (so that the early and late 
ISW can be neglected), it has been shown in Ref. Il62ll that for the s queezed lim it the effective /nl generated by second order 
gravitational effects on the CMB is /nl = —1/6 — 008(26*) (also see Ill50l4l52ll ). Here 6 is the angle between the short and the 
long modes. The angle dependent contribution comes from lensing. 



Effect of Cosmological parameter uncertainties 



Impact of uncertainties on the cosmol ogic al parameters effect the error bar on /nl- The effect of cosmological parameters 
have been discussed in Ref. ll25[|26ll28ill70ll . The cosmological parameters are determined using the 2-point statistics of the 
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CMB and therefor we expect the largest effect of /nl would come from those parameters which leave the CMB power spectrum 
unchanged while change the bispectrum. The expectation value of the estimator 

(f^L)-N 2^ Q,C,,C,3 ' ^^^^ 

changes with the change in cosmological parameters. Here Bi^i^i^ is the true CMB bispectrum. When changing the parameters 
the normalization should be c hanged to make the estimator unbiased. In general for a set of cosmological paramerets {p;}, 
the error in /nl is given by ll70ll 



5fNL = 

Here the average parameter values pi and their covariance matrix Cov{pi , pj ) can be determined using CMB-likelihood analysis 
tools. 

If the parameters are allowed to vary in the analysis then for WMAP this increases the Icr uncertainty in /nl by <5/nl° V/nl ~ 
16% for the local shape and 5/^1"*' / /nl ~ 14% for the equilateral shape. For Planck experiment the increases in Icr uncertainity 
is (^/nl^V/nl ~ 5% for local shape and (J/nl**' //nl ~ 4% for the equilateral shape. Most of the contribution to the error 
comes from three cosmological parameters, the amplitude of scalar perturbations A$, the tilt of the power spectrum of the saclar 
perturbations ns, and re-ionization optical depth t. 

For modes inside the horizon during reionization, the reionization optical depth r appears as a multiplicative factor e"'" in 
front of transfer function g^. For local model one of the mode is outside so the effect on bispectrum Ui^^^i^ ~ exp(— 2r)6^°^°^^ 

and for equilateral model all the modes are inside the horizon so ^Ij'^jls ~ ^^'P^^'^''')^7'ii2i3 ■ ''^^^ reduces to (5/n°l°' — — 2/nlt 
for local model and (J/nl"' — ^3/nlt for equilateral model. 

1 /2 

The effect of amplitude of perturbations can be seen by noting that the level of non-Gaussianity is given by /nl • . Hence 
the decrease (increase) in the amplitude of perturbations relax (tighten) the constraints on /nl- The effect of red tilt (n^ < 1) 
can be thought of as a reduction in power on at scales shorter than first peak and enhancement of power on larger scales. The 
effect of blue tilt is just opposite of red tilt. For local shape the limit on /nl becomes tighter proportional to A^^^^^ ||26|]. Note 
that Ref. il70ll show that the effect of cosmological parameters is negligible if the parameters are allowed to vary in the analysis 
and then marginalize over. 
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Cov(k,Pj) 



(58) 



Instrumental effects and distortions along the line of sight 

Here we point out that any cosmological or instrumental effect that can be modelled as a line of a sight CMB distortions of 
the primary CMB do not generate new bispectrum contribution. Although they can modify the the primordial bispectrum. A 
general model of line of sight distortions of the primary CMB are described in Ref. 1661,(681,071 

where the changes in the Stokes 

parameter of the CMB due to distortions along the line-of-sight can be written as 

6[Q ± iU]{n) - [a ± z2w](n)[Q ± iU]{h) + [h ± i/2](n)[Q T iU]{h) + [71 ± ^72](n)f (n) 

+ap(n) • V[Q ± iU]{h- a) + a[di ± id2]{-a)[di ± id2]f{h- a) + a^q{n)[di ± iSal^f (n; c^) + . . . (59) 

The first line captures the distortions in a single perfectly known direction n. The distortions in second line capture mixing 
of the polarization fields in a local region of length scale a around n. We Taylor expand the CMB fields Q, [/, and T around 
the point n and consider the leading order terms. Here Q, [/, and T stands for primordial (un-distorted) CMB fields. Since 
{QztiU) (n) is spin ±2 field, a(n) is a scalar field that describes modulation in the amplitude of the fields in a given direction n; 
a;(n) is also a scalar field that describes the rotation of the plane of polarization, (/i ± if 2) are spin ±4 fields that describe the 
coupling between two spin states (spin-flip), and (71 ± 172) (n) are spin ±2 fields that describe leakage from the temperature to 
polarization (monopole leakage hereon). Distortions in the second line of Eqn. ( |59] l. {pi ± P2), (^i ± ^2), and q are measured in 
the units of the length scale a. The field {pizLip2){n) is a spin ±1 field and describes the change in the photon direction; we will 
refer to it as a deflection field. Finally {di ± d2)(n) and q{h) describe leakage from temperature to polarization, {di ± d2){n) 
is spin ±1 field and we will refer to it as dipole leakage; q{h) is a scalar field that we will call quadrupole leakage. 

These distortions can be produced by various cosmological processes such as weak gravitational lensing of the CMB, screen- 
ing effects from patchy reionization, rotation of the plane of polarization due to magnetic fields or parity violating physics and 
various instrumental systematics such as gain fluctuations, pixel rotation, differential gain, pointing, differential ellipticity are 
also captured via line of sight distortions. All these distortions modify the primordial bispectrum as 
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(Pi' 

(2^ 



Q' Bit.M-i'M+nW^i^^ - nW^{l2 - I") + perm. 



(60) 



where is a window which depends on the type of distortion in consideration and tells how the primordial CMB bispectrum 
modes are coupled to the distortion field power spectrum Cf^. The effect of the distortions on the bispectru m is to sm ooth out 
the acousti c fea tures. These effects for the case of lensing have been shown to be small and can be neglected im [ml. 

In Ref. t l74|l the impact of the 1// noise and asymmetric beam on local has been found insignificant in the context of 

a Planck-like experiment. 



VI. OTHER PROBES OF NON-GAUSSIANITY IN THE CMB 



Although using the full bispectrum is the most sensitive cubic statistic other statistical methods may be sensitive to different 
aspects of non-Gaussianity and, more importantly, different methods have different systematic effects. Therefore it is important 
to study various probes. In this section we will discuss some of the methods which have been recently used or developed to test 
for primordial non-Gaussianities in the CMB. 



Trispectrum 

The four-point function in harmonic space is called trispectrum, which can be written as 

(ai,™,a,,™,a,3„3a,,„J = ( ^\ ) ( X ^ ^3 ) ^' ^^^^ 

where T^^^l'^ (L) is the angular averaged trispectrum, L is the length of a diagonal that forms triangles with li and I2 and with Z3 
and I4, and the matrix is the Wigner 3-j symbol. The trispectrum contains unconnected part, Tq, 

Tg\\\1{L) = (-l)'^+'V(2^i + 1)(2?3 + l)Ci,CiJi,iJi,ijLa 

+ {2L + l)Ci,Ci, [{-l)''+''+''5i,iJi,u + 5i,u5i2h] ■ (62) 

which comes from the Gaussian part of the perturbations, and the connected part Tc which contains non-Gaussian signatures. 
Using permutation symmetry, one may write the connected part of the trispectrum as 

Tc\±{L) = Pllll{L) + {2L + l)Y^ 

L' 

where 

Here, the matrix is the Wigner 6-j symbol, and i|^J^(i) is called the reduced trispectrum, which contains all the physical 
information about non-Gaussianities. For non-Gaussianity of local-type for which 

$(x) = $i(x) + /nl [$i - + 5NL$i , (65) 

both /nl and ^nl contribute to the trispectru m, bu t only /nl contributes to the bispectrum. Tripectrum based estimators for 
measuring f m, an d ^nl have been developed lll75l - ll79fl . For local template, the bispectrum nearly contains all the informa- 
tion on /nl 117711 . however if the non-Gaussianity is seen in bispectrum, trispectrum can s erve as a important cross-check. 
Generically for single field slow-roll models the trispectrum is small and un-observable ||I8C|] howeve r for more general single 
field models whenever the e quila teral bispectrum is large, the trispectrum is large as well lll8ll - ll83ll . For example for equi- 
lateral non-Gaussianity Ref. il84ll study how to tune the model parameters to get large trispectrum an d sm all bispectrum. For 
multi-field inflation one can construct models that predicts small /nl but large ^nl, for example Ref. Il85[l discusses the local 
form from a multi-field inflation, and briefly mentioned the condition in their class of models to get the large trispectrum and 
small bispectrum. Joint constraints on both /nl and gNL have the potential to add to the specificity of the search for primordial 
non-Gaussianity. For a giv en mo del these two numbers will often be predicted in terms of a single model parameter, such as a 
coupling constant, see e.g. II86II for the case of ek yprotic models. Using WMAP 5-year data, the constraints on ^nl using the 
trispectrum are -7.4 < .gNL/10^ < 8.2 at 2cr il87ll . 
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Minkowski Functionals 



Minkows ki Functio nals (MFs) describe morphological properties (such as area, circumference, Euler characteristic) of fluctu- 
ating fields il88l4l9l]l . For a d-dimensional fluctuating field, /, the A:-th Minkowski Functionals of weakly non-Gaussian fields 

in, vlf\v) for a given threshold v = //cto can be written as lll92l[T93ll 



ao + OK)^, (66) 



where o-q = if^)^^^ is the variance of the fluctuating field, Hn{v) are the Hermite polynomials, LOk = 7r'""/^/r(fc/2 + 1), and 
finally S*''^ are the "skewness parameters" defined as 

EE (67) 

c(2) _ 3d ((V/) ■ (V/)(VV)) 

^ = ^w^) ^ ' ^ ^ 

which characterize the skewness of fluctuating fields and their derivatives. Here ai characterizes the variance of the fluctuating 
field and is given by 

^l^^Y.^'2i+m{i+^)rcj (70) 

For CMB, for which d = 2 and / = AT/T, the skewness parameters are il94ll 

c(l) _ 3 hih + 1) + hjh + 1) + ^3(^3 + 1) pmim2m3^mim2m3T;rA M/ W/ 



where Bi^i^i^ is the CMB bispectrum, Wi represents a smoothing kernel which depends on the experiment beam and 5™;™;^™^ 
is the usual Gaunt function. 

Since MFs can be determined as weighted sum of the bispectrum, they contain less information than the bispectrum. MFs can 
still be useful because they perhaps suffer from different systematics, though they are less specific to primordial non-Gaussianity 
since they measure a smaller number of independent bispectrum modes. Also, the bispectrum is defined in Fourier (or harmonic) 
space while the MFs are defined in real space. Limits on non-Gaussianity of local-type from the MFs of the WMAP 5-year 
temperatur e data are —70 < /nl < 91(2(t) I195i1 . The MFs from the Planck temperature data should be sensitive to /nl ~ 20 
at 1(7 level ll96ll in contrast to bispectrum which is sensitive to /nl ^ 5 at Icr level. Note that polarization data further improves 
the sensitivity. 



Wavelets 

Several studies have used wavelet representations of the WMAP maps to search for a non-Gaussian signal il97l4203]l . In most 
of these studies, wavelets were used as a tool for blind searches of non-Gaussian anomalies in a basis with resolution in both scale 
and location. However, in some more recent studies, wavelets were tuned to look for non-Gaussianity of a particular type. In the 
context of searches for primordial non-Gaussianity of local type, wavelet based estimators for /nl have been built by extracting a 
signature of local non-Gaussianity that is cubic in the wavelet coefficients from simulations of non-Gaussian skies and searching 
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for this signature in data. This ability to calibrate on a set of simulations makes the wavelet approach very flexible. While not 
optimal in a least-squared sense, using a wavelet representation can be thought of as a generalized cubic statistic with a different 
weighting scheme to the optimal bispectrum estimator Using such estimators therefore provides a useful exploration of nearly 
optimal cubic estimators similar to the full bispectrum estimator Any believable detection of non-Gaussianity should be robust 
to such changes in the analysis. Similarly, contaminating non-Gaussianity from astrophysical and instrumental systematics will 
propagate through the analysis in a different way to the bispectrum-based analysis. 

There are several constraints on local /nl using wavelet based estimators. For example, using the COBE data the constraints 
are |/nl| < 2200(lcr) i204| |. Using an estimator based on the skewness of the wavelet coefficients, Mukherjee and Wang 
constrain the /nl value for WMAP 1-yr data obtaining /nl = 50 ± 160(2(t) ll29ll . Using an extension of the pr eviou s estimator 
by combining wavelet coefficients at different contiguous scales, Curto et al. obtain — 8 < /nl < lll(2tT) l205ll . Rec ently, 
using a generalized third order estimator based on the wavelet coefficients, Curto et al. obtain — 18 < /nl < 80(2(7) l206ll . 



Needlet Bispectrum 

Needlets are a family of spherical wavelets which are localized and asymptotically unco rrected |207ll208ll . The needlet 
based statistics as been considered for testing Gaussianity and isotropy (for example see Ref. 12091 - 1215 1. Using the bispectrum 
of needlet coefficient, the constraints on non-Gaussianity of local-type using WMAP 5-year data yields /nl = 73 ± 62(2(7) 
121611211^ . As is clear, the needlet based bispectrum is not as sensitive as the CMB bispectrum discussed in Sec. lIVI however 
again in the event of detection the needlet based methods can be calibrated on simulations and represent a different weighting 
scheme for handling the sky mask and anisotropic noise. Finally, needlets and wavelets allow for the possibility to analyze 
spatially localized regions in the sky. 



Probing non-Gaussianity using Bayesian statistics 



A somewhat different approach to searching for non-Gaussianity is provided by the Bayesian approach. Here, the starting 
point is an explicit physical or statistical model for the data and the goal is to evaluate the posterior density of the parameters of 
the model and/or the relative probability of the Gaussianity and non-Gaussianity. 

On large scales, in the Sachs- Wolfe regime, one can simplify the Bayesian approach by modeling directly the temperature 
anisotropy. Rocha et al. (2001) l218ll discuss a Bayesian exploration of a model where each spherical harmonic coefficient 
is drawn from a non-Gaussian distribution. In this regime, the simple form of th e non -Gaussian potential for the local model 
Eq. ( fTSb also translates into a simple model for the temperature anisotropy. Ref. l219ll develop several results for it, including 
an analytical expression of the evidence ratio of the Gaussian and non-Gaussian models. At the level of current data, this 
approximation is too restrictive, since most of the information about /nl is contained near the resolution limit of the experiment, 
where most of the measured perturbation modes are concentrated. 

A full implementation of a physical non-Gaussian model must include the effect of Boltzmann transport. In the context of 
local non-Gaussianity, the model equation Eq. ( fTSl ) suggests that a full Bayesian treatment may be feasible. At the time of 
writing, no fully Bayesian analysis for local /nl has been published. The effort has focused on developing approximations to 
the full Bayesian problem. 

Using a perturb ative analysis, Ref. l22Clll relate the frequentist bispectrum estimator to moments of the Bayesian posterior 
distribution. Ref. l22lll described approximations to the full Bayesian treatment that simplify the analysis for high signal-to- 
noise maps and compare these to the full Bayesian treatment for a simple 1-D toy model of non-Gaussian sky where this analysis 
is feasible. 



VII. SUMMARY 



The physics of the early universe responsible for generating the seed perturbations in the CMB is not understood. Inflation 
which is perhaps the the most promising paradigm for generating seed perturbations allow for vast number of inflationary models 
that are compatible with data based on 2-point statistics like CMB power spectrum. Moreover, the alternatives to inflation such 
as cyclic models are also compatible with the data. Characterizing the non-Gaussianity in the primordial perturbations has 
emerged as probe for discriminating between different models of the early universe. Models based on slowly rolling single field 
produce undetectable amount of non-Gaussianity. Single field models without the slow roll can generate large (detectable with 
future experiments) non-Gaussianities but 1) can not produce large non-Gaussianity of local type unless inflation started with 
with excited vacuum state 2) if non-Gaussianity is produced it would naturally be as bispectrum while higher order such as 
trispectrum can be generated, it requires fine tuning. 
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TABLE II: Summary of constraints on local non-Gaussianity 



/^fr'±2aeiTor 



Year 



data 



Method 



|/nl| < 1500 
I/nlI < 1900 
39.5 ± 97.5 
/nl < 5400 
47± 74 

30 ±84 
32 ±68 
87 ±62 

110 ± 910 
10.5 ±80.5 
51 ±60 

'"-950 

58 ±46 
38 ±42 

31 ±49 
-315 ± 705 

11 ±47.4 

32 ±42 



2002 
2003 
2003 
2004 
2005 
2006 
2006 
2007 
2007 
2008 
2008 
2008 
2009 
2009 
2009 
2009 
2009 
2010 



COBE 
MAXIMA 
WMAP 1-year 

VSA 
WMAP 1-year 
WMAP 3-year 
WMAP 3-year 
WMAP 3-year 

Boomerang 
WMAP 3-year 
WMAP 5-year 
ARCHEOPS 
WMAP 3-year 
WMAP 5-year 
WMAP 5-year 
BOOMERanG 
WMAP 5-year 
WMAP 7-year 



Bispectrum sub-optimal 
Bispectrum sub-optimal 
Bispectrum sub-optimal 
Bispectrum sub-optimal 
Bispectrum sub-optimal-vl 
Bispectrum sub-optimal 
Bispectrum sub-optimal-vl 
Bispectrum near-optimal 
Minkowski Functionals 
Minkowski Functionals 
Bispectrum near-optimal 
Minkowski Functionals 
Bispectrum optimal 
Bispectrum optimal 
Spherical Mexican hat wavelet 
Minkowski Functionals 
Skewness power spectrum 
Bispectrum optimal 



Komatsu et al. [2221 
Santos et al. [2231 
E. Komatsu et al. [231 
Smith et al. [2241 
Creminelli et al. [251 
Spergel et al. [241 
Creminelli et al. [261 
Yadav and Wandelt [2^ 
De Troia et al. [2251 
C. Hikage et al. [1951 
Komatsu et al. [511 
Curto et al. 2008 [2261 
Smith et al. [I3I1 
Smith et al. [1311 
Curto, A et. al. [2061 
R Natoli et al. [2271 
Smidt, Joseph et al. [2281 
Komatsu et al. [1321 



The bispectrum of the CMB is one of the most promising tool for connecting the non-Gaussianities in the cosmic microwave 
background and the models of inflation. Bispectrum-based estimator which saturates Cramer-Rao bound has been developed 
and well characterised using non-Gaussian Monte-Carlos. Other statistics although not as sensitive to non-Gaussianity as an op- 
timally weighted bispectrum estimator, do provide independent checks and have different systematics. While Bayesian analysis 
has been applied in the context of non-Gaussianity analysis, this still appears to be an open area for fruitful investigation. 

Given the importance of detecting primordial non-Gussianity, it is crucial to characterise any non-primordial sources of non- 
Gaussianities. We describe several sources of non-Gaussianities such as from second order anisotropies after last scattering 
surface and during recombination. 

With Planck launched and taking data, we look forward to the next few years as an exciting time in the exploration of 
primordial non-Gaussianity in the cosmic microwave background. 
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